Asymmetry and decoherence in a double-layer persistent-current qubit 
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Superconducting circuits fabricated using the widely used shadow evaporation technique can con- 
tain unintended junctions which change their quantum dynamics. We discuss a superconducting flux 
qubit design that exploits the symmetries of a circuit to protect the qubit from unwanted coupling 
to the noisy environment, in which the unintended junctions can spoil the quantum coherence. We 
present a theoretical model based on a recently developed circuit theory for superconducting qubits 
and calculate relaxation and decoherence times that can be compared with existing experiments. 
Furthermore, the coupling of the qubit to a circuit resonance (plasmon mode) is explained in terms 
of the asymmetry of the circuit. Finally, possibilities for prolonging the relaxation and decoherence 
times of the studied superconducting qubit are proposed on the basis of the obtained results. 



I. INTRODUCTION 



Superconducting (SC) circuits in the regime where the 
Josephson energy Ej dominates the charging energy Ec 
represent one of the currently studied candidates for a 
solid-state qubit 0. Several experiments have demon- 
strated the quantum coherent behavior of a SC flux 
qubit 0, IE 3 , and recently, coherent free- induction de- 
cay (Ramsey fringe) oscillations have been observed 
The coherence time extracted from these data was 
reported to be around 20 ns, somewhat shorter than ex- 
pected from theoretical estimates 0, El El- I R morc 
recent experiments [ljj, it was found that the decoher- 
ence time T2 can be increased up to approximately 120 ns 
by applying a large dc bias current (about 80% of the 
SQUID junctions' critical current). 

A number of decoherence mechanisms can be impor- 
tant, being both intrinsic to the Josephson junctions, e.g., 
oxide barrier defects \\\\ or vortex motion, and exter- 
nal, e.g., current fluctuations from the external control 
circuits, e.g., current sources 0,0, H E3 El Here, we 
concentrate on the latter effect, i.e., current fluctuations, 
and use a recently developed circuit theory 01 to ana- 
lyze the circuit studied in the experiment |5| ■ 

The SC circuit studied in Ref. @ (see Fig. UJ is de- 
signed to be immune to current fluctuations from the 
current bias line due to its symmetry properties; at zero 
dc bias, Ib = 0, and independent of the applied magnetic 
field, a small fluctuating current 61 b (t) caused by the fi- 
nite impedance of the external control circuit (the current 
source) is divided equally into the two arms of the SQUID 
loop and no net current flows through the three-junction 
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FIG. 1: Schematic of the circuit. Crosses denote Josephson 
junctions. The outer loop with two junctions L and R forms 
a dc SQUID that is used to read out the qubit. The state of 
the qubit is determined by the orientation of the circulating 
current in the small loop, comprising the junctions 1, 2, and 
3, one of which has a slightly smaller critical current than 
the others. A bias current Ib can be applied as indicated for 
read- out. 



qubit line. Thus, in the ideal circuit, Fig. ^ the qubit is 
protected from decoherence due to current fluctuations 
in the bias current line. This result also follows from a 
systematic analysis of the circuit However, asymme- 
tries in the SQUID loop may spoil the protection of the 
qubit from decoherence. The breaking of the SQUID's 
symmetry has other very interesting consequences, no- 
tably the possibility to couple the qubit to an external 
harmonic oscillator (plasmon mode) and thus to entan- 
gle the qubit with another degree of freedom |13| . For an 
inductively coupled SQUID a small geometrical 

asymmetry, i.e., a small imbalance of self-inductances in 
a SQUID loop combined with the same imbalance for the 
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FIG. 2: Schematics of Josephson junctions produced by 
the shadow evaporation technique, connecting the upper with 
the lower aluminum layer. Shaded regions represent the alu- 
minum oxide. 



mutual inductance to the qubit, is not sufficient to cause 
decoherence at zero bias current [T^ . A junction asym- 
metry, i.e., a difference in critical currents in the SQUID 
junctions L and R, would in principle suffice to cause 
decoherence at zero bias current. However, in practice, 
the SQUID junctions are typically large in area and thus 
their critical currents are rather well-controlled (in the 
system studied in Ref. |10| . the junction asymmetry is 
< 5%), therefore the latter effect turns out to be too 
small to explain the experimental findings. 

An important insight in the understanding of deco- 
herence in the circuit design proposed in is that it 
contains another asymmetry, caused by its double layer 
structure. The double layer structure is an artifact of 
the fabrication method used to produce SC circuits with 
aluminum/ aluminum oxide Josephson junctions, the so- 
called shadow evaporation technique. Junctions pro- 
duced with this technique will always connect the top 
layer with the bottom layer, see Fig. [3 Thus, while cir- 
cuits like Fig. ^ can be produced with this technique, 
strictly speaking, loops will always contain an even num- 
ber of junctions. In order to analyze the implications of 
the double layer structure for the circuit in Fig. ^ we 
draw the circuit again, see Fig. OJa), but this time with 
separate upper and lower layers. Every piece of the upper 
layer will be connected with the underlying piece of the 
lower layer via an "unintentional" Josephson junction. 
However, these extra junctions typically have large areas 
and therefore large critical currents; thus, their Joseph- 
son energy can often be neglected. Since we are only 
interested in the lowest-order effect of the double layer 
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FIG. 4: External circuit attached to the qubit (Fig. that 
allows the application of a bias current Ib for qubit read- 
out. The inductance L s h and capacitance C s h form the shell 
circuit, and Z(uj) is the total impedance of the current source 
(Ib)- The case where a voltage source is used to generate a 
current can be reduced to this using Norton's theorem. 



structure, we neglect all unintentional junctions in this 
sense; therefore, we arrive at the circuit, Fig.[3fb), with- 
out extra junctions. We notice however, that this result- 
ing circuit is distinct from the 'ideal' circuit Fig. ^ that 
does not reflect the double-layer structure. In the real 
circuit, Fig. |3fb), the symmetry between the two arms 
of the dc SQUID is broken, and thus it can be expected 
that bias current fluctuations cause decoherence of the 
qubit at zero dc bias current, Ib = 0. This effect is par- 
ticularly important in the circuit discussed in |5lll(| since 
the coupling between the qubit and the SQUID is dom- 
inated by the kinetic inductance of the shared line, and 
so is strongly asymmetric, rather than by the geometric 
mutual inductance Q which is symmetric. Our analysis 
below will show this quantitatively and will allow us to 
compare our theoretical predictions with the experimen- 
tal data for the decoherence times as a function of the 
bias current. Furthermore, we will theoretically explain 
the coupling of the qubit to a plasma mode in the read- 
out circuit (SQUID plus external circut), see Fig. at 
Ib = 0; this coupling is absent for a symmetric circuit. 

This article is organized as follows. In Sec. |Hj we de- 
rive the Hamiltonian of the qubit, taking into account 
its double-layer structure. We use this Hamiltonian to 
calculate the relaxation and decoherence times as a func- 
tion of the applied bias current (Sec. IIIIII and to derive 
an effective Hamiltonian for the coupling of the qubit to 
a plasmon mode in the read-out circuit (Sec. II V|) . Fi- 
nally, Sec. IVl contains a short discussion of our result and 
possible lessons for future SC qubit designs. 



FIG. 3: (a) Double layer structure. Dashed blue lines repre- 
sent the lower, solid red lines the upper SC layer, and crosses 
indicate Josephson junctions. The thick crosses are the in- 
tended junctions, while the thin crosses are the unintended 
distributed junctions due to the double-layer structure, (b) 
Simplest circuit model of the double layer structure. The 
symmetry between the upper and lower arms of the SQUID 
has been broken by the qubit line comprising three junctions. 
Thick black lines denote pieces of the SC in which the upper 
and lower layer are connected by large area junctions. 



II. HAMILTONIAN 

In order to model the decoherence of the qubit, we 
need to find its Hamiltonian and its coupling to the en- 
vironment. The Hamiltonian of the circuit Fig. 03[b) can 
be found using the circuit theory developed in Ref. |12| . 
To this end, we first draw the circuit graph (Fig. |SJ) and 
find a tree of the circuit graph containing all capacitors 



FIG. 5: The network graph of the circuit, Figs, and 2] 
Dots indicate the nodes, lines the branches of the graph; an 
arrow indicates the orientation of a branch. Thick lines la- 
beled Ji denote an RSJ element, i.e., a Josephson junction 
shunted by a capacitor and a resistor. Lines labeled Li and 
Ki denote inductances, Z ext the external impedance, includ- 
ing the shell circuit of Fig. 2] and Ib is the current source. 



and as few inductors as possible (Fig. EJ). A tree of a 
graph is a subgraph containing all of its nodes but no 
loops. By identifying the fundamental loops ^3 m the 
circuit graph (Fig. |SJ we obtain the loop submatrices 
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The chord (L) and tree (K) inductance matrices are 
taken to be 



_ / L/2 M/4 \ 
^ _ ^ M/4 L'/2 J ' 



■'K 



/ L/2 M/4 M t 
M/4 L'/2 
V Mi Li 



(3) 

where L, L' , and Li are, respectively, the self-inductances 
of the qubit loop in the upper layer, the SQUID, and 
qubit loop in the lower layer, and M and Mi are the 
mutual inductances between the qubit and the SQUID 
and between the upper and lower layers in the qubit loop. 
The tree-chord mutual inductance matrix is taken to be 



M/4 
M/4 



(4) 



The Hamiltonian in terms of the SC phase differences 
ip = ((px,cp2,(p3,(pL,(pR) across the Josephson junctions 
and their conjugate variables, the capacitor charges Qc, 
is found to be Hal 
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FIG. 6: A tree of the circuit graph, Fig. A tree is a 
subgraph connecting all nodes, containing no loops. Here, 
the tree was chosen to contain all capacitors Ct (from the RSJ 
elements) and as few inductors Ki as possible (see Ref. ^3)- 
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where the Josephson inductances are given by Lj-i = 
$o/27r/ c; i, and L c -i is the critical current of the i-th junc- 
tion. In Eq. ©, we have also introduced the effective 
self-inductances of the qubit and SQUID and the effec- 
tive qubit-SQUID mutual inductance, 
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(7) 
(8) 
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and the coupling constants between the bias current and 
the qubit and the left and right SQUID phases, 



kT^I + L'/L + 2M/i)(l - 2MJL), (10) 



1 1 

2 ~ 26' 

with the definitions 



m R 
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25' 



(11) 



k = l+AL t (L + L' + 2M)/L 2 + 2(L' + M ~2Mi)/L 
-(M + 2Mi) 2 /L 2 , (12) 
5 = k/(1 + M/L + 2M l /L)(l-2M i /L). (13) 

The sum ip± + ip2 + (f3 is the total phase difference 
across the qubit line containing junctions J\, J2, and J3, 
whereas ipL + <fiR is the sum of the phase differences in the 
SQUID loop. Furthermore, C = diag(C, C, C, C , C) is 
the capacitance matrix, C and C being the capacitances 
of the qubit and SQUID junctions, respectively. 

The working point is given by the triple (/,/', Lb), 
i.e., by the bias current Lb, and by the dimensionless 
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external magnetic fluxes threading the qubit and SQUID 
loops, / = 27r$ 2 ./$ and /' = 27r$^/$ - We will work 
in a region of parameter space where the potential U(tp) 
has a double-well shape, which will be used to encode the 
logical qubit states 0) and 1). 

Ine classical equations of motion, including dissipa- 
tion, are 
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where convolution is defined as (/ * g)(t) — f f(t — 
T)g(r)dT. The vector m is given by 
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m = A(mQ,m Q ,m Q ,m L ,m R ), 



(15) 



and A is chosen such that |m| = 1. For the coupling 
constant /i, we find 

/i = k- 2 L- 4 (Z(L + L' + 2M) 2 (L + M -2M t ) 2 

+ (2 Li (L + L' + 2M)+L (i'-2M j ) (16) 
—M (M + 2 Mi)) 2 + (L 2 + 2 Li (L + L' + 2 M) 

+L (L' + 2M-2 Mi) - 2 M t (M + 2 A/,)) 2 ) . 

The kernel K in the dissipative term is determined by 
the total external impedance; in the frequency domain, 



Z{w) 



with the impedance 

Z(u>) = Z cxt (uj) + iujL int , 
where we have defined the internal inductance, 



(17) 
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-L(2L' 2 + 2L'M - M 2 + AM Mi - %M 2 )) , 



and where 



Z B (u) 



iujC sh 



iloL= 



(20) 



is the impedance of the external circuit attached to the 
qubit, including the shell circuit, see Figs. 14151 For the 
parameter regime we are interested in, L m t ~ 20 pH, lu < 
10 GHz, and Z > 50 ft therefore LoL lnt ^ |^cxt|, an d we 
can use Z(u>) ~ Z ex t(w). 

We numerically find the double-well minima tp and 
(p 1 for a range of bias currents between and 4 /iA and 
external flux f /2tt between 1.33 and 1.35 and a qubit 
flux around f/2n ~ 0.5 (the ratio ///' = 0.395 is fixed 
by the areas of the SQUID and qubit loops in the cir- 
cuit). The states localized at (p Q and (p x are encod- 
ing the logical |0) and |1) states of the qubit. This 



FIG. 7: Decoupling (red solid) and symmetric (blue dashed) 
curves in the (Ib, /') plane, where Ib is the applied bias cur- 
rent and /' = 27T<E>^./$o is the dimensionless externally ap- 
plied magnetic flux threading the SQUID loop. Both curves 
are obtained from the numerical minimization of the potential 
Eq. ©. The decoupling line is determined using the condi- 
tion m ■ A(p = 0, whereas the symmetric line follows from the 
condiction e = 0. 



allows us to find the set of parameters for which the 
double well is symmetric, e = U(<p Q ) — E/(Vi) = 0. 
The curve /*(Ib) on which the double well is symmet- 
ric is plotted in Fig. Qualitatively, /*(/b) agrees 
well with the experimentally measured symmetry line 
[Tof. but it underestimates the magnitude of the vari- 
ation in flux /' as a function of Ib- The value of Ib 
where the symmetric and the decoupling lines intersect 
coincides with the maximum of the symmetric line, as 
can be understood from the following argument. Tak- 
ing the total derivative with respect to Ib of the relation 
e = U( V o;r(I B ),I B )- U{fp x ;r{I B \I B ) = on the 
symmetric line, and using that y x are extremal points 
of U, we obtain n • A<pdf*/dI B + (27r/$ )m ■ A<p = 
for some constant vector n. Therefore, m • Aip = (de- 
coupling line) and n • Aip ^ implies df* jdls = 0. 

For the numerical calculations throughout this pa per , 
we use the estimated experimental parameters from [l(l 
E3, L = 25 pH, V = 45 pH, M = 7.5 pH, L t = 10 pH, 
Mi = 4pH, I C , L = I C . R = 4.2 [iA, and I c . x = I c , 2 /a = 
I 3 = 0.5 /iA with a ~ 0.8. 



III. DECOHERENCE 

The dissipative quantum dynamics of the qubit will be 
described using a Caldeira-Leggett model [Tj| which is 
consistent with the classical dissipative equation of mo- 
tion, Eq. JT1J. We then quantize the combined system 
and bath Hamiltonian and use the master equation for 
the superconducting phases ip of the qubit and SQUID in 
the Born-Markov approximation to obtain the relaxation 
and decoherence times of the qubit. 



5 



A. Relaxation time Ti 

The relaxation time of the qubit in the semiclassical 
approximation |15| is given by 



. ( 21) 

where Ay> = tp — tp 1 is the vector joining the two minima 
in configuration space and 



E = VA 2 + e 2 



(22) 



is the energy splitting between the two (lowest) eigen- 
states of the double well and A is the tunnel coupling 
between the two minima. We will evaluate T\ on the 
symmetric line where e = 0, and therefore, E — A. At 
the points in parameter space (Ib, /') where m- Acp van- 
ishes, the system will be decoupled from the environment 
(in lowest order perturbation theory), and thus T± — > oo. 
From our numerical determination of <p and tp 1 , the de- 
coupling flux /', at which m • A(p — 0, is obtained as a 
function of Ib (Fig. CJ. From this analysis, we can infer 
the parameters (Ib, /') &t which T\ will be maximal and 
the relaxation time away from the divergence. In prac- 
tice, the divergence will be cut off by other effects which 
lie beyond the scope of this theory. However, we can fit 
the peak value of T\ from recent experiments |ldj | with 
a residual impedance of R rcs ~ 3.5 Mil which lies in a 
different part of the circuit than Z (FigEJ. We do not 
need to further specify the position of i? rcs in the circuit; 
we only make use of the fact that it gives rise to an ad- 
ditional contribution to the relaxation rate of the form 
Eq. H21(l but with a vector m rcs ^ m, with m res • Acp ^ 
on the decoupling line. Without loss of generality, we 
can adjust R rcs such that m res ■ Aip = 1. Such a residual 
coupling may, e.g., originate from the subgap resistances 
of the junctions. The relaxation time Ti obtained from 
Eq. (|21|l as a function of Is along the symmetric line 
e = (Fig. [7J with a cut-off of the divergence by R les is 
plotted in Fig. El along with the experimental data from 
sample A in In Fig. [SJ we also plot T\ (theory and 
experiment) as a function of the applied magnetic flux 
around the symmetric point at zero bias current. For the 
plots of T\ in Figs. |H1 and El we have used the experimen- 
tal parameters A/h = 5.9 GHz, Z(E) ~ Z cxt (E) = 60 Q, 
and T = 100 mK. 



B. Decoherence time T-2 

The decoherence time T2 is related to the relaxation 
time T\ via 

1 1 1 



T 2 T, 



<t> 



2Ti' 



(23) 



where denotes the (pure) dephasing time. On the 
symmetric line /' = J*(Ib) (see Fig. CJ, the contribu- 
tion to the dephasing rate TT 1 of order Rq/Z vanishes, 




FIG. 8: Theoretical relaxation time Ti (solid line) as a func- 
tion of the applied bias current Ib, along the symmetric line 
(Fig. 0. The value of I B where Ti diverges coincides with 
the intersection of the symmetric line with the decoupling 
line in Fig. the divergence is removed in the theory curve 
by including a residual impedance of R las = 3.5 MQ. The ex- 
perimentally obtained data for sample A in Ref.El are shown 
as triangle symbols with error bars. 



where Rq = e 2 /h w 25.8 Ml denotes the quantum of 
resistance. However, there is a second-order contribu- 
tion ex (Rq/Z) 2 , which we can estimate as follows. The 
asymmetry e = U((f ) — of the double well as a 

function of the bias current Ib at fixed external flux /' 
can be written in terms of a Taylor series around I B , 



e(I B ) = eo + ei&Zu + e 2 5I B + 0(SI E 



(24) 



where SlB(t) = Ib(P) ~I% IS the variation away from the 
dc bias current I B . The coefficients ej(/g) can be ob- 




1.345 



FIG. 9: Theoretical relaxation time Ti (solid line) as a 
function of the applied magnetic flux /' = $^./$o at zero 
bias current, Ib = 0, around the symmetric point, e = 0. 
Experimentally obtained data for sample A in Ref. [l(J are 
shown as triangle symbols with error bars. The theory curve 
from the semiclassical Ti formula, Eq. 121H . is expected to 
be valid in the range |e| < A, which corresponds roughly to 
1.33 < /' /2tt < 1.34. Experimental points outside the plot- 
ted range of /' where the theory curve is not expected to be 
valid, are not shown. 
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The calculated relaxation and decoherence times T\ 
and T 2 agree well with the experimental data p"(i| in their 
most important feature, the peak at Ib ~ 2.8 fj,A. This 
theoretical result does not involve fitting with any free 
parameters, since it follows exclusively from the inde- 
pendently known values for the circuit inductances and 
critical currents. Moreover, we obtain good quantitative 
agreement between theory and experiment for T\ away 
from the divergence. The shape of the T\ and T 2 curves 
can be understood qualitatively from the theory. 

IV. COUPLING TO THE PLASMON MODE 



FIG. 10: Theoretical relaxation, pure dephasing, and deco- 
herence times Ti, T^, and T2 as a function of applied bias 
current Ib, along the symmetric line (Fig. |TJ| . As in Fig. [HI 
we have included decoherence from a residual impedance of 
R ICB = 3.5 Mfi. 



tained numerically from the minimization of the poten- 
tial U, Eq. ©. The approximate two- level Hamiltonian 
<tz in its eigenbasis is then, up to 0(5I B ), 



H 



A 

2 



4A 



eoei 
~2A~' 



Ua 



£q£2 
2A 



<r z 5I 2 



(25) 
l» (26) 



where A = A + £q/2A. On the symmetric line, eo = 
0, the term linear in SIb vanishes. However, there is a 
non- vanishing second-order term cx e\ that contributes to 
dephasing on the symmetric line. Without making use of 
the correlators for 51 B , we know that the pure dephasing 
rate T^ 1 will be proportional to 6i(Ib) 4 which allows us 
to predict the dependence of on Ib- A discussion of 
the second-order dephasing within the spin-boson model 
can be found in |l6j| . However, in order to explain the 
order of magnitude of the experimental result 10] for 
correctly, the strong coupling to the plasma mode may 
also play an important role [lfllr7T |. The result presented 
here cannot be used to predict the absolute magnitude of 
Tqj, but we can obtain an estimate for the dependence of 
on the bias current Ib via €\{Ib) = de/dls obtained 
numerically from our circuit theory, via 



T^(Ib) « T7\0) 



^i(Ib) 
ei(0) 



(27) 



where for dimensional reasons we can write the propor- 
tionality constant in terms of a zero-frequency resistance 
Rq and an energy Q (note, however, that this corresponds 
to one free parameter in the theory), IT (0)/ei(0) ~ 
2u> 3 /R%A 2 . For the plots of T and T 2 in Fig. M we 
have used the resistance Ro — 1450 f2 and have chosen 
u>/2tt w 1 THz to approximately fit the width of the T 2 
curve. The relaxation, dephasing, and decoherence times 
Ti, Ta, and T 2 are plotted as a function of the bias cur- 
rent I B in Fig. and Fig.lTUl 



In addition to decoherence, the coupling to the exter- 
nal circuit (Fig. can also lead to resonances in the 
microwave spectrum of the system that originate from 
the coupling between the qubit to a LC resonator formed 
by the SQUID, the inductance L s h and capacitance C s h 
of the "shell" circuit (plasmon mode) . We have studied 
this couplin g q uantitatively in the framework of the cir- 
cuit theory [12] , by replacing the circuit elements Ib and 
Z in the circuit graph by the elements L s h and C s h in 
series, obtaining the graph matrices 



'CL 



(- 1 1 

-1 1 

-1 1 

-1 

-1 



°\ 




1 





-1 1 
= I 0-11), (28) 
-1 1 



V 0-1/ 



where the last row in Fcl corresponds to the tree branch 
C s h and the rightmost column in both Fc l and F^- j, cor- 
respond to the loop closed by the chord L s h- Neglecting 
decoherence, the total Hamiltonian can be written as 



H = Hs + Ush + Hs,sh, 



(29) 



where Hs, defined in Eq. (JSJ), describes the qubit and 
SQUID system. The Hamiltonian of the plasmon mode 




FIG. 11: Plasma frequency cj s h as a function of the applied 
bias current Ib ■ The variation is due to the change the effec- 
tive in Josephson inductances as Ib is varied. 
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can be brought into the second quantized form 

**-^ + (s)'£-**(" + 5)- <»> 

by introducing the resonance frequency w s h = l/V-^tCsh, 
the total inductance (where the SQUID junctions have 
been linearized at the operating point) L t ~ L s h + L'/4 + 
L'j / (cos(ip l) + cos(ipn)), and the creation and annihila- 
tion operators W and b, via 

<t>o V zGsh^sh V Kq 



with the impedance = \jL t j C s h . For the coupling 
between the qubit/SQUID system (the phases <p) and the 
plasmon mode (the phase </? s h associated with the charge 
on C sh , <3 sh = C sh $o^sh/2vr, we obtain 

^ sh= (£) i^" 1 -^ (32) 

where m is given in Eq. JTHJ and M sh w L sh + Z//4 (the 
exact expression for M s h is a rational function of L s h and 
the circuit inductances which we will not display here). 
Using Eq. (|3*T|) and the semiclassical approximation 

m • ip rj — icr z m • Ac^> + const., (33) 

we arrive at 




Ib [M] 

FIG. 12: Coupling constant A between the qubit and the 
plasmon mode. The coupling disappears at the crossing with 
the decoupling line (Fig. [TJ , i.e., when m • A(p = 0. 

qubit and the plasmon mode. In particular, recently ob- 
served side resonances with the sum and difference fre- 
quencies E ± oj s h |13| can be explained in terms of the 
coupled dynamics, Eq. (|36[) . Also, it should be possible 
to tune in-situ the coupling to the plasmon mode A at 
will, using pulsed bias currents. 



V. DISCUSSION 



W s ,sh = Aa 2 (6 + 6t), 
with the coupling strength 




m ■ A(p. 



(34) 



(35) 



Note that this coupling vanishes along the decoupling line 
(Fig. and also rapidly with the increase of L s h- 

The complete two-level Hamiltonian then has the well- 
known Jaynes-Cummings form, 



H = Aa x + ea z + hu sh 



+ \a z (b + tf). (36) 



For the parameters in Ref. fToL C s h = 12 pF and L s h = 
170 pH, we find uj sh » 2tt x 2.9 GHz (see Fig. Ell and 
Z sh = 5 H, thus ^Z sh /R Q fa 0.01. Note that the de- 
pendence of the Josephson inductance (and thus of L t 
and w s h) on the state of the qubit leads to an ac Stark 
shift term cx a z Wb which was neglected in the coupling 
Hamiltonian Eq. 1)36(1 . 

We find a coupling constant of A w 210 MHz at Ib = 0. 
The coupling constant as a function of the bias current 
Ib is plotted in Fig. El The relatively high values of A 
should allow the study of the coupled dynamics of the 



We have found that the double-layer structure of SC 
circuits fabricated using the shadow evaporation tech- 
nique can drastically change the quantum dynamics of 
the circuit due to the presence of unintended junctions. 
In particular, the double-layer structure breaks the sym- 
metry of the Delft qubit [j| (see Fig. , and leads to re- 
laxation and decoherence. We explain theoretically the 
observed compensation of the asymmetry at high Ib [To| 
and calculate the relaxation and decoherence times T\ 
and T 2 of the qubit, plotted in Fig. El We find good 
quantitative agreement between theory and experiment 
in the value of the decoupling current Ib where the re- 
laxation and decoherence times T\ and T 2 reach their 
maximum. In future qubit designs, the asymmetry can 
be avoided by adding a fourth junction in series with the 
three qubit junctions. It has already been demonstrated 
that this leads to a shift of the maxima of T\ and T 2 close 
to Ib = 0, as theoretically expected, and to an increase 
of the maximal values of T\ and T 2 [T(| • 

The asymmetry of the circuit also gives rise to an inter- 
esting coupling between the qubit and an LC resonance in 
the external circuit (plasmon mode) , which has been ob- 
served experimentally [13| , and which we have explained 
theoretically. The coupling could potentially lead to in- 
teresting effects, e.g., Rabi oscillations or entanglement 
between the qubit and the plasmon mode. 
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